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Abstract. We will give an introduction to the concept of quarkyonic matter, presenting an 
overview of what is meant by this term in the literature. We will then argue that the quarkyonic 
phase, as defined in the original paper, is a percolation-type phase transition whose phase 
transition line is strongly curved in ps — N c space, where N c is the number of colors and ps 
the baryon density. With a toy model estimate, we show that it might be possible to obtain a 
percolating but confined phase at N c — 3, Nf — 2 at densities larger than one baryon per one 
baryon size. We conclude by discussing how this phase can be observed at FAIR. 



A concept recently receiving quite a lot of theoretical attention is that of "quarkyonic matter" . 
Introduced in [lj, it can be best defined as matter where temperature is below deconfinement, 
T < T c , but quark chemical potential is of the order of the confinement scale (inverse of the 
baryon size), hq > Aqcd> resulting in a density of a baryon per baryonic volume or higher. 
This regime is not easily accessible through controlled theoretical expansions: Effective hadronic 
models, since they expand in powers of Q/Aqcd> where Q is the momentum exchange, can be 
reasonably expected to diverge. Lattice calculations are affected by the sign problem, which 
generally can only be controlled at \xq <C T. 

Hence, expectations for matter in this regime are very model-specific. Popular ideas included 
a confined yet chirally restored phase [3] , chirally inhomogeneus phases [U [6] , a deconfined 
but chiral-broken plasma of constituent quarks [7] and a large-iV c generalization of nuclear 
matter [8j [10]. Ref. [Tj has generically argued that in this regime confinement coexists 
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Figure 1. Diagrams at arbitrary N c responsible for gluon antiscreening and quark- hole 
screening. 



with pQCD and asymptotic freedom. The essence of this argument is rather simple and 
qualitative: Assuming that at finite chemical potential deconfinement appears when the quark- 
hole screening overpowers the gluon-gluon antiscreening, we can obtain how the chemical 
potential required for deconfinement varies with N c . It turns out that at leading order (Fig. [TJ, 
^confinement ^ / ]V^VjAqcd ' ^ higher order (n loops) the N c dependence is parametrically 
unchanged. Thus, looking at momentum space, we expect hq ~ Aqcd matter to be firmly 
confined. 



In configuration space, on the other hand, we have a baryon with iV c quarks per baryonic 

— 1/3 

size. Hence, quarks of neighboring baryons are distant by N c — > 0. Thus, at least at large N c , 
asymptotic freedom should somehow coexist with confinement. A natural picture of this, found 
in [T], is an asymptotically free gas of quarks moving in a mean field of baryons. Alternatively, 
the bulk of the Fermi surface is asymptotically free but its excitations are confined. 
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Figure 2. The influence of percolation of quark wavefunctions inside a dense distribution of 
hadrons (shown here as potential wells). 

How relevant is this reasoning to our N c = 3, Nf = 2 world? As we have already known 
for a few decades, the N c — > oo expansion is qualitatively remarkably good for a lot of QCD 
phenomenology: That mesons are weakly-interacting quasi-particles, baryons look like semi- 
classical "skyrmion" states of mesonic fields, and planar diagrams dominate are all features of 
the large N c expansion. On the other hand, some features of QCD thermodynamics are not 
even qualitatively well reproduced: At large N c , deconfinement is a first order phase transition 
rather than a cross-over. At large chemical potential, binding energy of nuclear matter is of the 
order of the baryon mass ~ _/V c Aqcd (rather than ~ 10~ 3 -/V c Aqcd *C m^, as in our world). 

Thus, the existence of phase transitions in N c , as well as T and fiQ, is virtually assured. Such 
transitions can not obviously be studied through experimental data, but nevertheless remain 
pretty interesting theoretically: The large-iV c technique is one that has been extensively used 
for phenomenology as well as theory. The identification of the large--/V c limit with a classical 
dual gauge theory makes this topic more interesting still, since the phase transition in N c 
is a quantum-gravitational one in its dual description. Unfortunately, as confinement shows, 
determining that a phase transition exists from a fundamental theory can be a very subtle issue. 

In the rest of this work, we shall use percolation [T2] as a "toy model" to investigate 
a transition between normal nuclear matter and the quarkyonic phase which bears a lot of 
similarity with that described in pp. Percolation is a simple yet surprisingly deep and universal 
paradigm: It requires an infinite network of "links" , each with a probability, invariant in space, 
to be "on" or "off". The crucial question is the typical length of a connected cluster of "on" 
links as a function of the individual link probability^ Somewhat subtly, this link's dependence 
on probability is not monotonic: At the critical "percolation" probability (typically not very 
high, O (1) /Neighbours), this length diverges. Considering this as a correlation length, it is 
clear percolation defines a second order phase transition, an insight which is confirmed by the 
appearance of critical behavior at the percolation threshold. The abstractness and generality of 
percolation makes it very appealing as a toy model, since any local propagator in an extended 
medium should give rise to percolation. Specifically, by analogy with the well-known conductor- 
insulator transition in QED, percolation signals the regime where one can not consider quark 
wavefunctions to be effectively localized to a baryon since quantum tunneling will have the 
potential to make quarks propagate through arbitrary distances at equilibrium (Fig. [2]). 

To investigate the curvature in ps-N c space, we posutulate a generic interaction propagator, 
respecting what we know about Yang-Mills theory: Due to confinement, it has to fall off 
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Figure 3. The dependence of the percolation critical number of colors N c on A, at a density 
of one baryon per baryon size and hexagonal packing [TT] . 
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Figure 4. Two-dimensional plots in the (baryon density,AT c ) plane of the percolation curves 
(lines without data points) vs. the deconfmement curves at various temperatures (8 = T/T c ), 
for Nj = 1,3 (resp. left and right plots). 



at distances higher than Aq CD in configuration space, and the interaction probability has to 

~A/iV c : 




in units of Aqcd> it is natural that ry ~ O (1). Since at that scale A becomes strong, it is natural 
that that it is A ~ O (1) too. As long as these requirements are satisfied, the shape of f(x) 
(0-function in position space, Gaussian, etc) is irrelevant. Fig. [3] |TTj shows the dependence of 
the critical number of colors -/V c for a regular hexagonal packing on A, at fixed baryon density. 
It is clear that at the chemical potential of (Iq ~ Aqcd> indicated in [lj, iV c = 3 is not enough 
to trigger a percolation transition, which requires N c ~ O(10). However, the possibility of 

quarkyonic matter at 1 < hq /Aqcd < \J N c /Nf remains open if the effective number of flavors 
Nf < 3. Note that, since the mass of the strange quark ~ Aqcd> it is far from clear whether the 
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Figure 5. Spectral function of dileptons in a hadron gas, in a pQCD, and in quarkyonic matter. 



above hyerarchy can be maintained, since it is uncertain whether the strange flavor is "light" 
for our purposes. 

To further investigate whether the quarkyonic phase can appear, we varied density as well 
as N c [12] . Note that, while this is an obvious extension, technically it is somewhat difficult 
because the definition of "nearest neighbor" is somewhat sloppy with varying density (the nearest 
neighbor at one density can be farther than the next-to-nearest neighbor at the other). To obtain 
the critical percolation probability, we integrated out physical baryons into a square lattice, and 
performed a renormalization-group analysis. The details of this procedure are available in |12j . 

Unlike percolation, baryonic density necessary for deconfinementgoes up with N c because 
of the hq ~ yj 1 N c /Nf Aqcd requirement. Hence, the "true high-A^ limit" is one where 
deconfinement happens at higher [iq than percolation. In this limit, there will be both a 
confined percolating phase and a deconfined phase. At the low N c limit, percolation would 
require a higher \jlq than deconfinement. Since percolation requires baryonic wavefunctions to 
be present, this means percolation newer occurs. Where is our world in this hyerarchy? In Fig. H] 
|12j we plotted, as a function of N c , /j,q,T, the density of an ideal baryon gas at deconfinement. 
As can be seen, a sliver of percolating but confined phase diagram remains. While the "ideal" 
approximation for baryon gas density is rough in this limit, this calculation shows percolation 
is possible in our world. Hence, an experimental search for it might be necessary. 




Figure 6. Influence of the baryonic dynamics on the decomposed quark- hole spectral function 
in quarkyonic matter. 



What would be an experimental signature for the production of "quarkyonic matter"? This 
question was not really answered by pQ; however, the coexistence, in quarkyonic matter, 



of pQCD degrees of freedom with baryonic mean fields makes exploration of the spectral 
function by dilepton detection a promising avenue to investigate. Consider a reaction such 
as q + hole — > 7* — > l + l~ whose rate V will be given by 

r(M) = J M ((pi - p 2 ) 2 ) f(pi)f(p2)d 3 Pl d 3 p 2 5 ((p? - p£) 2 - M 2 ) ; (2) 

the assumption that quark and hole degrees of freedom be quark-like means A4 (...) should not 
be too different from the quark-quark scattering matrix elements calculated within pQCD. 

The distribution functions f(p± 2) will not however be those of a free gas but will contain a 
mean field reflecting the distribution of baryons in the system. Equivalently, the Eigenstates of 
the unperturbed quarkyonic matter will not be free particles (i.e., the spectral function will not 
be flat in M), but rather Eigenstates of the Hamiltonian 

(baryons) 

H = k 2 + m 2 q + V mi (x ) , V mf (x,y)= J2 $ mf (x - (3) 

i 

where the mean field baryonic potential $ m f will contain dips centered around the position of the 
baryons in configuration space (see Fig. [2]). If baryons were immobile and at regular intervals, 
as we assumed in Fig. [3] and we would recover Bloch's theorem for QED conductors: the 

1 /3 

spectral function of dileptons would experience a dip at M ~ p^ , as shown in Fig. [5l due 

1/3 

to the disappearance of quark wavefunctions having a momentum ~ p^, . Such dips would be 
distinctive for quarkyonic matter, since a hadron gas would give peaks for the production of 
resonances and a perturbative QGP would give a relatively flat spectral function (Fig. [5]). 

Density of inhomogeneities and the movement of baryons are likely to at least somewhat spoil 
this picture. However, the dynamics of baryons is calculable by semiclassical dynamical models 
such as uRQMD [13], and shows regularities such as transverse and longitudinal expansion. Also, 
experimentally we can access not just the invariant mass of the spectral function but its full 
vector decomposition My jT ^ in rapidity-transverse-azimuthal space. Since in quarkyonic matter 
quarks are not point particles but non-trivial Eigenfunctions of the distribution of baryons, 
the shape of the M y ^ r ^ distribution is likely going to be non-trivial, and reflecting macroscopic 
dynamics. To calculate it, one would need to solve the Dirac equation for dips centered around 
a classical distribution of baryons (whose positions are given by RQMD), and then fold the 
solutions with pQCD interaction matrices. This work is currently in progress, and we expect to 
get quantitative predictions for dilepton distributions from quarkyonic matter. 

In conclusion, we have used a percolation-based model to argue for the existence of a dense 
"quarkyonic" phase in which quark wavefunctions can percolate to arbitrary distances, yet 
confinement persists and baryons are well-defined objects. We have shown that this phase is 
almost certain to exist in the limit N c 3> Nf, and is not excluded in our N c = 3,Nf = 2 
world. We have proposed dilepton spectral functions as a possible experimental verification 
of this phase, and await developments at low-to-intermediate energy heavy ion collisions, such 
as the FAIR program. G. T. acknowledges the financial support received from the Helmholtz 
International Center for FAIR within the framework of the LOEWE program (Landesoffensive 
zur Entwicklung Wissenschaftlich-Okonomischer Exzellenz) launched by the State of Hesse. 
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